Ultracompact Generation of Continuous- Variable Cluster States 
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We propose an experimental scheme that has the potential for large-scale realization of continuous- 
variable (CV) cluster states for universal quantum computation. We do this by mapping CV cluster- 
state graphs onto two-mode squeezing graphs, which can be engineered into a single optical para- 
metric oscillator (OPO). The desired CV cluster state is produced directly from a joint squeezing 
operation on the vacuum using a multi-frequency pump beam. This method has potential for ultra- 
compact experimental implementation. As an illustration, we detail an experimental proposal for 
creating a four-mode square CV cluster state with a single OPO. 
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One-way quantum computation [l[ is a promising form 
of quantum computing (QC). Unitary gates are imple- 
mented by performing single-qubit measurements on a 
highly entangled "cluster state" Q • Such states are par- 
ticular cases of graph states Q and are characterized by 
this ability to coherently manipulate quatum informa- 
tion using only local measurements and classical feed- 
forward of measurement results. One often discusses 
QC in terms of two-dimensional systems (qubits) but 
continuous quantum variables can also be used 0, Hi- 
Continuous- variable (CV) cluster states @, 0] generalize 
the concept of cluster states to quantum systems with 
continuous degrees of freedom by use of the quadratures, 
amplitude q = (a + a') and phase p = —i(a — at), of a 
boson field Q. The preparation procedure exactly mir- 
rors that for qubit cluster states, using the correspon- 
dence between the Pauli and continuous Weyl-Heisenberg 
groups and can be described in an analogous way 
to the graph state formalism Q for qubit cluster states: 
First, prepare each mode (represented by a graph vertex) 
in a phase-squeezed state, approximating a zero-phase 
eigenstate (analog of Pauli-AT eigenstates). Then, ap- 
ply a quantum nondemolition (QND) interaction of the 
form exp(iqjqk) (analog of controlled-phase Cz) to each 
pair of modes (j, k) linked by an edge in the graph. All 
QND gates commute (as do Cz gates), so the full multi- 
mode QND operator to be applied is exp(±q T Aq), where 
q = (qi, . . . ,otv) T is a vector of amplitude quadrature 
operators, and A is the adjacency matrix for the graph. 
The resulting cluster state satisfies, in the limit of infinite 
squeezing, the relation Q 



p - Aq — > 



(1) 



as a simultaneous zero-eigenstate of the N components 
of p — Aq. Being infinitely squeezed, this state is un- 
physical (it has infinite energy) but we can approximate 
its Wigner function by any Gaussian strongly squeezed 
along the indicated quadratures. Defining a "CV cluster 



state" as any of these approximating Gaussians is iden- 
tical to the definition used in Ref. [9( and more general 
than that studied in Refs. 0, 0]. While only Gaussian 
operations are needed to create CV cluster states from 
the vacuum, using them for universal QC requires that at 
least one single-mode non-Gaussian measurement (such 
as photon-number resolving detection) be available Q ■ 

Although convenient theoretically, the above proce- 
dure is not optimal for experimental implementation be- 
cause the QND gates contain in-line squeezers. It can be 
spectacularly simplified 9] by use of the Bloch-Messiah 
reduction [lfj , which transforms any Gaussian operation 
into the canonical form of a set of single-mode squeezers 
(e.g., optical parametric oscillators — OPOs) sandwiched 
between two multimode interferometers. With a vac- 
uum input, the initial interferometer is irrelevant and any 
Gaussian CV N-mode cluster can be formed by N single- 
mode vacuum squeezers (easier to implement than in-line 
ones) followed by a network of 0(N 2 ) beam splitters, i.e. 
a quadratically large, stable interferometer [9J. Recently, 
a four- mode linear cluster state was demonstrated fllj |. 

In this Letter, we show that it is, in fact, possible 
to integrate all single-mode OPOs into one multimode 
OPO, pumped by an 0(N 2 )-mode field, and to eliminate 
the beam splitter network completely. This is equally 
resource-efficient as the proposal in Ref. 0] but the com- 
plexity has been shifted from a stabilized 0(A^ 2 )-element 
interferometer (unwieldy for large N) to the nonlinear 
medium of a single OPO and the frequency content of 
the pump beam. This scheme is much more compact 
and it has interesting prospects for scalability because 
it effectively represents the quantum entangled version 
of an optical frequency comb: as is well known, a fem- 
tosecond laser effectively compactifies ~ 10 5 phase-locked 
continuous- wave lasers into a single beam T^, HI- One 
of us showed that such a comb can be transformed into a 
GHZ state / dx \x) 1 ■ ■ ■ \x) N , where the subscripts denote 
consecutive comb lines, using a complete network of con- 
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current nonlinear interactions 14J, [l5j, LL6| , and the nonlin- 
ear medium required to create four-mode entanglement 
in a single OPO has already been demonstrated lTf. En- 
gineering concurrent nonlinear interactions between an 
arbitrary number of modes is a complicated problem but 
it is now solvable in the general case by use of general- 
ized quasi-phase-matching in photonic quasicrystals [l8| . 
This enables arbitrarily difficult nonlinear interactions 
(e.g., simultaneous generation of the second, third, and 
fourth optical harmonics, all in different directions [HI]) 
to be engineered in a single OPO. 

The central result of this paper is a mathematical 
connection between CV cluster states and two-mode 
squeezing (TMS) graph states flfjj ]. CV cluster-state 
graphs have vertices representing phase-squeezed states 
and edges corresponding to QND unitary interactions 
and are the ones we wish to implement for one-way 
QC 0- TMS graphs have vertices representing vacuum 
inputs and edges denoting individual terms in the multi- 
mode squeezing Hamiltonian 
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where G denotes the adjacency matrix of the graph. We 
prove that the two are related: any CV cluster state with 
a bipartite graph can be created by applying a single 
multimode squeezing Hamiltonian of the form of Eq. ([2|) 
and any such Hamiltonian generates some CV cluster 
state. We detail how to create a square cluster using this 
method with current technology. 

Given a target CV cluster state, our goal is to effect 
a transformation on the quadrature operators such that 
Eq. (TTJ) holds for the new quadratures. We first collect q 
and p into a column vector x = (<?i, . . . , qN,Pi, ■ ■ ■ ,Pn) T ■ 
Gaussian transformations on the vacuum in Hilbert space 
correspond to symplectic linear transformations on this 
vector in the Heisenberg picture [l9|. We denote by U a 
the symplectic transformation corresponding to a unitary 
that creates a CV cluster state from the vacuum. The 
level of overall squeezing is represented by a > 0, which 
should be as large as possible. From Eq. |T]) we have 



{-A l)U a x ^0, 



(3) 



where the block matrix above is N x 2N, U a is 27V x 2N, 
and Xq is the vector of quadrature operators representing 
the vacuum state. The arrow denotes the limit a — > oo. 

We have some additional freedom in Eq. ([3]). After 
the transformation U a is applied, we can perform arbi- 
trary phase shifts for each individual mode at the out- 
put, which we will represent with the matrix T. This 
is a passive transformation on the state, which can be 
effected simply by reinterpreting the output modes (i.e., 
no change to the physical apparatus used to create the 
state is required). Therefore, we have that 



is sufficient to conclude that U a can be used to create 
a CV cluster state with adjacency matrix A from the 
vacuum. 

As we will now show, if A represents a bipartite graph, 
we can always do this with a multimode squeezing Hamil- 
tonian. By definition, the nodes of a bipartite graph are 
partitioned into two sets such that all graph edges link 
one set to the other. These graphs are also known as 
two- colorable graphs because the two sets (and the nodes 
each contains) can be assigned different colors. Bipartite 
graphs include the square lattice graph of arbitrary size, 
which is universal for QC, and any of its subgraphs. Star 
graphs (of any size) are also bipartite, with the node at 
the center being one color and the rest a different color. 
As a counterexample, the triangle graph (and, more gen- 
erally, any graph with an odd cycle in it) is not bipartite. 

Consider a multimode squeezing Hamiltonian given by 
Eq. J2]), where G is the (as yet, undefined) adjacency 
matrix for a TMS graph. Writing U a as the Heisenberg 
matrix corresponding to exp(— ictH) gives 
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A large (but finite) value of a is required for a useful CV 
cluster state. Although previous work 14, 15, [l6| has 
emphasized uniformly weighted TMS graphs with no self- 
loops, at this point the only restriction we are going to 
place on G is that it be symmetric and full-rank. Experi- 
mental requirements will favor some G's over others but, 
since any G is in principle possible to implement [3] , we 
will not impose any additional restrictions at this point. 

With these requirements we can write G as the differ- 
ence of two positive semidefinite matrices that are mutu- 
ally orthogonal. By this we mean G = G+ — G_, where 
G± > and G±G T = G T G± = 0. We write G° ± for 
the Moore-Penrose pseudoinverse of G± , which (for sym- 
metric matrices) is obtained by inverting all the nonzero 
eigenvalues of G±. Then, G" 1 = G% - G°_. The pro- 
jectors onto the positive and negative subspaces of G are 
P± = G±G° ± = G° ± G±. Recalling Eq. ©, we need both 
the positive and negative exponentials of G in the limit 
of large a. In the positive (negative) case, such an oper- 
ation will magnify all the positive (negative) eigenvalues 
of G and zero out all of G's negative (positive) eigenval- 
ues. To write this concisely, we start with the fact that 
e~ aG± — > P T for large a, since all of the nonzero eigen- 
values of G± get sent to zero (since G± > 0) while the 
zero eigenvalues get raised to 1. This gives 

e ±aG = e -"G Te aG± ^ p ±e "G± = G±G ^ e "G± _ (g) 

By suitably numbering nodes, the adjacency matrix for 
any bipartite graph can be written as 
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A = 
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where A is L x (JV — L). Instead of using colors, we 
will label the first L modes by + and the rest by — be- 
cause the number of each will correspond to the number 
of positive and negative eigenvalues of G, respectively. 
Recalling Eq. (TJ|, we will use the phase-shift freedom in 
T to rotate all of the — modes by —tt/2 and leave the 
others unchanged. This gives 
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where I± is the identity matrix on the ± modes and zero 
on =F modes and the identity blocks and zero blocks are 
sized appropriately, according to the dimensions of Aq. 
Plugging Eqs. ©, ©, and © into Eq. g]) gives the 
following sufficient condition for cluster state creation: 



I A \ fG+G° + e aG + 
-Al JO oil 





G-G°_e aG - 
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Keeping in mind that the first matrix is N x 2N, while 
the second is 2N x 2N, this condition is fulfilled if 



(-AS I)G + = Q 



and 



(I Ac) G_ = . (10) 



These requirements are satisfied by choosing 



B (I A ) , G_ 
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where B, C > are arbitrary symmetric positive defi- 
nite matrices. This also illustrates our earlier point that 
labeling the sets of nodes as + and — reflected their con- 
nection to the number of eigenvalues of G having each 
sign. Thus, a CV cluster state with a bipartite adja- 
cency matrix A satisfying Eq. ([7]) can be created with a 
TMS Hamiltonian of the form of Eq. with 



G = 
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For a given A, this is the most general G that satisfies 
Eq. ([9]), since B and C encompass all possible rotations 
of the eigenvectors and scalings of the eigenvalues that 
preserve the partitioning defined by Eq. Ijl0|> . With A 
fixed, the freedom to choose the G that is easiest to im- 
plement experimentally is found solely within the choices 
of B and C. 

This is not the most general solution to the overar- 
ching problem, however. There is no reason a priori 
that we should have a completely fixed A for a given 
CV cluster state that we wish to create. While all QND 



interactions in the original formulation Q of CV clus- 
ter states for QC had the same strength, this is not 
necessary. A weighted adjacency matrix A corresponds 
to variable-strength QND interactions for the edges of 
the graph. This introduces squeezing and/or reversal 
(q —>■ aq, p — > p/a, where a is the edge weight) to the 
Gaussian correction term that accumulates after each 
measurement. While very low (or very high) weights 
would lead to difficulty resolving the quantum state af- 
ter being heavily squeezed, for weights ~ ±1, both the- 
oretically and practically speaking, all of the quantum 
information is still preserved under single-mode measure- 
ments made on the cluster. Allowing A to be weighted 
gives additional degrees of freedom to the problem, al- 
lowing us even greater freedom in optimizing the experi- 
mental viability of the multimode squeezing Hamiltonian 
used to make our cluster state. 

A corollary to this result is that any multimode squeez- 
ing Hamiltonian of the form of Eq. ([2|) that has a full-rank 
G generates some weighted bipartite CV cluster state (af- 
ter appropriate single- mode phase shifts). To see this, 
write G in terms of its eigendecomposition G = VvV T , 
where v is a diagonal matrix of eigenvalues and V is an 
orthogonal matrix. Using elementary column operations, 
up to a possible renumbering of the output modes, we can 
always transform V into the form of the first matrix in 
Eq. (fT2"| . The target form always exists because it is the 
simultaneous column-reduced echelon form for the pos- 
itive and negative subspaces of G and G is assumed to 
be full-rank. These column operations, since they act 
separately on the two subspaces, can be represented by 
an invertible block-diagonal matrix M acting from the 
right, such that VM = ( Jr f°)- The transpose of this 

matrix, M T , acting from the left, represents the same ac- 
tion as row operations on V T . With M being invertible 
and block-diagonal, we can choose a B, C > such that 
M~ x vM~ T = (-q Thus, we can always write 



G = VM{M- 1 vM- T )M T V 1 
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for some particular Aq . Comparing this with Eq. (|12j) , we 
can immediately extract Aq and use Eq. (|7|) to write A in 
terms of it. This completes the proof. We therefore also 
know that any multimode squeezing Hamiltonian gen- 
erates a weighted bipartite CV cluster state (generally 
with a different graph A) as long as the TMS adjacency 
matrix G is full-rank. 

Intuitively, what's happening with this correspondence 
is that H from Eq. ([2]) is used to squeeze the vacuum 
along N joint quadratures (since G is full-rank) with 
overall squeezing strength a. In general, these states 
are not CV cluster states because they do not satisfy 
Eq. |1| for any choice of A in the large-a limit. What we 
have shown is that by partitioning the resulting output 
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modes into two groups (corresponding to the number of ± 
eigenvalues of G) and phase-shifting one of those groups 
by — 7r/2, we can always transform the output from the 
multimode squeezer into a CV cluster state, satisfying 
Eq. (fTJ) for some choice of A as a becomes large. Our 
derivation requires that A be bipartite for this to work. 

As an example, let G be the complete graph on four 
nodes. This generates a GHZ state 14| whose quadrature 
operators satisfy 91+92 + 93 + 94 — > 0, Pi — P2 —* 0, 
P1—P3 — * 0, andpi— P4 — > (and any linear combinations 
thereof). Phase-shifting mode 1 (although any mode will 
do) by — 7r/2 means that now — pi + 92 + 93 + 94 — > 0, 
91 — P2 —* 0, 91 —p3 — > 0, and 91—^4 — > 0, which satisfies 
Eq. (fT]) with A being the star graph on four nodes, with 
node 1 in the center. This property generalizes: G being 
the complete graph on N nodes creates an TV-mode GHZ 
state, which is equivalent to an TV-mode star-graph CV 
cluster state after phase-shifting one of the output modes 
by — 7r/2. The shifted mode becomes the central node in 
the star. This mimics the case for qubits ]3j, although 
the analogy is not exact since G and A represent different 
types of graphs (TMS and CV-cluster, respectively). 

Star graphs are not universal for QC, however. We'd 
like to achieve a procedure for generating a square-lattice 
(or other QC-universal) CV graph with a single OPO. 
Such a graph is bipartite, so a corresponding G can be 
constructed to create it and, in principle, can be quasi- 
phase- matched in a single photonic quasicrystal [181 ]. A 
significant step in this direction is the creation of a CV 
cluster state with a square graph from a single four-mode 
OPO: one can, indeed, show the remarkable result, 
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FIG. 1: Experimental implementation of a square CV cluster 
state using a single OPO based on a periodically poled bire- 
fringent crystal, such as KT1OPO4 (see also Ref. \L7\)- Left: 
the cluster graph (A) after a phase-shift of modes 3 and 4 by 
— it/2; dashed line denotes a negative weight; all magnitudes 
are l/y2. Right: the experimental proposal. Top arrows are 
OPO modes, bottom arrows are pumps, and all have polar- 
ization directions y or z along the crystal axes. Nonlinear 
interactions simultaneously phasematch (first letter is pump) 
yyz (open circles), yzy (filled circles), zzz (open squares), 
and zyy (filled squares). The OPO cavity resonance con- 
ditions and crystal birefringence ensure that no other OPO 
mode can be coupled to these four modes. Note the crucial 
importance of the 7r-shifted pump, —z. 

cluster states using only one OPO and no beam-splitter 
network. 
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Notice that A is weighted so that one of the edges (sides of 
the square) has an opposite interaction sign to the three 
others and all have magnitude l/y/2. A (nonunique) gen- 
erating G is identical and immediately implementable us- 
ing current technology, in fact, using the existing nonlin- 
ear crystal [l7| designed to produce the four-party CV 
GHZ state (Fig. [I]). Defining a as in Eq. (JSJ), the vari- 
ance in each of the components of p — Aq for this state 
is 2e~ 2a units of vacuum noise. Since these vanish as 
a — > 00, this is a valid square-graph CV cluster state. 

In conclusion, we have shown that any continuous- 
variable cluster state with a bipartite graph can be gen- 
erated from the application of a single multimode squeez- 
ing Hamiltonian. We also have shown that all multimode 
squeezing Hamiltonians that have a full-rank two-mode 
squeezing adjacency matrix correspond to a weighted bi- 
partite continuous-variable cluster state, generally with 
a different graph. While as resource-efficient as the most 
efficient scheme currently known Q, these results arc 
important for experiments because they provide a pow- 
erfully scalable means of generating continuous- variable 



[1] 
[2] 
[3] 

[4] 

[5] 
[6] 

[7] 

[8] 

[9] 

[10] 

[11] 
[12] 
[13] 
[14] 



Electronic address: nmen@princeton.edu 

Electronic address: opfister@virginia.edu 

R. Raussendorf and H. J. Briegel, Phys. Rev. Lett. 86, 

5188 (2001). 

H. J. Briegel and R. Raussendorf, Phys. Rev. Lett. 86, 
910 (2001). 

M. Hein et al., Proceedings of the International School of 
Physics "Enrico Fermi", on Quantum Computers, Algo- 
rithms and Chaos, 162 (2006), quant-ph/0602096. 
S. Lloyd and S. L. Braunstein, Phys. Rev. Lett. 82, 1784 
(1999). 

S. D. Bartlett et al, Phys. Rev. Lett. 88, 097904 (2002). 



L. Braunstein, Phys. Rev. A 73, 032318 
Phys. Rev. Lett. 97, 110501 



et al. 



J. Zhang and S 
(2006). 

N. C. Menicucci 
(2006). 

D. F. Walls and G. 
(Springer, Berlin, 1994). 
P. van Loock, C. Weedbrook 
ph/0610119. 

S. L. Braunstein, Phys. Rev. A 71, 055801 (2005) 
X. Su et al, Phys. Rev. Lett. 98, 070502 (2007). 
T. W. Hansen, Rev. Mod. Phys. 78, 1297 (2006). 
J. L. Hall, Rev. Mod. Phys. 78, 1279 (2006). 
O. Pfister et al, Phys. Rev. A 70, 020302 (2004). 



J. Milburn, Quantum Optics 



and M. Gu, quant- 



■5 



[15] A. S. Bradley et ai, Phys. Rev. A 72, 053805 (2005). 

[16] O. Pfister, quant-ph/0701104. 

[17] R. C. Pooser and O. Pfister, Opt. Lett. 30, 2635 (2005). 

[18] R. Lifshitz, A. Arie, and A. Bahabad, Phys. Rev. Lett. 



95, 133901 (2005). 
[19] R. Simon, N. Mukunda, and B. Dutta, Phys. Rev. A 49, 
1567 (1994). 



